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0^ ' Abstract. We show that strongly localized wave lunctions occur around classical bifurca- 

0^ , tions. Near a saddle node bifurcation the scaling of the inverse participation ratio on Planck's 

constant and the dependence on the parameter is governed by an Airy function. Analytical 
estimates are supported by numerical calculations for the quantum kicked rotor. 
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Many mechanisms of localization in dynamical systems have been discussed. The 
saturation in the spreading of a wave packet in the quantum kicked rotor can be 
' connected to Anderson localization ||l|, ||. Other effects, such as trapping due to 

cantori Q or in elliptic islands, have a classical origin. The scarring phenomenon, 
an increase in intensity near classical periodic orbits, seems to be intermediate [|], ^. 
, Scars are most prominent on orbits that are weakly unstable or marginally stable, 

' as in the stadium billiard. This suggests that near a classical bifurcation, where the 

, semiclassical amplitudes diverge, they should be most prominent. As we will show here 

this is indeed the case: the classical slowing down in phase space near a bifurcation 
causes an enhanced eigenfunction amplitude, but there is also a quantum interference 
between two orbits that brings in quantum modulations. Specifically, we here study 
the localization of the wave functions in the quantized standard map near saddle node 
Q ' bifurcations of period one orbits. 

, The quantum kicked rotor is described by a finite unitary matrix which in the 

momentum plane wave representation is given by H 
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■ The dimension N of the momentum basis is related to Planck's constant by ?i = An/N. 

The parameters k and t in (^ are connected to h and the classical kicking strength 
parameter, K,hy k = K/t and r = Att/N, where the latter choice ensures the absence 
of dynamical localization. In the classical limit the quantum kicked rotor corresponds 
to the standard map defined on a torus, 

p„+i = p„ + ifsin i9„ mod 47r, i3„+i = + p„+i mod 27r . (2) 

The classical standard map depends on K only. For K > the phase space is mixed 
and at K ^ 0.93996... the last KAM torus breaks up. For K > 5 the system is 
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essentially fully mixing with exponentially small stable islands. However, even for suf- 
ficiently large K values the classical phase space is very complex. Of particular interest 
here are parameter values near a saddle node bifurcation where pairs of periodic orbits 
are created or destroyed. Such bifurcations carry large semiclassical weights since the 
second derivative of the action vanishes and the semiclassical amplitude diverges. The 
effect of such bifurcations on the spectral statistics has been much studied [||. Here 
we focus on the wave functions. 

The simplest examples of saddle node bifurcations arise when K passes through 
K* = Att. The two pairs of stable and unstable period one orbits are created. Fig. |l| 
shows the phase space structure for a if-value slightly above the bifurcation, K/A-k = 
1.0216. 
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Fig. 1 Classical phase space near the bifurcation for K/4n — 1.0216. The bifurcation takes 
place at {&o = 7i"/2, po — 27r) 

The quantity we use to detect localization in wave functions is the average inverse 
participation ratio of the wave functions, i.e. the fourth moment of their components 
Cai = (^Iq^) of the a's eigenstate in some orthonormal basis \i) summed over all eigen- 
states and normalized by the value expected within random matrix theory for the 
orthogonal ensemble |l2|, 

a i 

Thus, if all wave functions are ergodically spread over all basis states, L = 1. Localized 
states will show an L > 1. This is a global quantity that is related to the return 
probability which can be calculated in the following way. The transition probability 
from state i to state j after m kicks is Pij (m) = |Tr(Aij[/'")p, where Aij is a projector 
Aij = from the i-th to the j-th basis state. The localization measure L is 

nothing but the average long time return probability, L = l/3X]i^Mi where Pa — 
liniM^oo X]m=i Piii.^^) the return probability to state 

The quantity L can also be estimated using the semiclassical procedures proposed 
in 1^ for matrix elements. Specifically, let An = be the projector on the i-th 

basis state. Then = \ ~ 'fa) the matrix element weighted density 

of states. The semiclassical expression for has two terms, pA^ = PoiAu) + Posc, 
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where the first term is the sniooth phase space average and comes from the Thomas- 
Fermi density of states, whereas the latter part can be expressed semiclassically as 
a sum over periodic orbits. The weight of each periodic orbit p also depends on the 
integral of the observable along the points of the orbit. The fourth moment that 




Fig. 2 Measure of localization in the vicinity of the saddle-node bifurcation at K/4n = 1 
for different values of = 201, 401 and 801. The curves have been rescaled according to the 
uniformization expected from Eq. 

enters the expression for L is a measure of the fluctuations of the matrix elements 
(a|Aii|a) = |(z|q;)P — \cai\'^ around their average and can hence be calculated from 
the square of the sum over periodic orbits. Within the diagonal approximation we 
therefore obtain 

P 

where the Wp contain the stability information about the orbit and the the con- 
tribution from the observable When all orbits are uniformly hyperbolic the 
random matrix value is recovered. However, in our situation near a bifurcation we 
get diverging amplitudes Wp which, when expressed in a uniform approximation, are 
replaced by an Airy function. For the fixed points near K* /Att =integer, we can write 
Pose = Yl'p ~^tbif, where the sum extends over all orbits except the ones near the bi- 
furcation. The weights tp are given in a conventional form |^ whereas the bifurcating 
orbit contribution |l§ is t^f « Ai„fN^^'^Ai[N'^/^{K-K*)]. Therefore the fiuctuations 
are expected to show a contribution which is a square of an Airy function 

L-l^N^/^Ai^[bN^/^{K~K*)]. (5) 

In Fig. H we present the variation of L as a function of K around the saddle-node 
bifurcation of period-one orbits at K* /Air — 1 for different values of N. For K < K* 
the random matrix behavior is recovered with increasing N . Beyond K* Airy-type 
oscillations appear. Both the K and the L values have been rescaled according to the 
A^-scaling in (|^). After rescaling the first peaks for different TV the curves fall on top 
of each other. Higher than second moments of the eigenvector components show an 
even enhanced effect provided the distribution function of |ccip develops a tail due to 
the presence of localized states. 
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As an illustration in Fig. |^ we show the state with the largest inverse participation 
ratio at the value of K/Att = 1.0216 for N = 401. This state is localized at p = 2-k. In 
^^-representation the state is indeed localized around d — tt/2 and d — 37r/2, right at 
the pair of orbits that is born in the saddle-node bifurcation. Due to the uncertainity 
principle its spatial extent is larger than in p-representation. 
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Fig. 3 A localized state at K = 1.0216 for A'^ = 401 in p~ (left panel) and i9— representation 
(right panel). The constant line indicates the mean value expected from random matrix 
theory. 

In summary, we have analyzed the localization of quantum eigenstates in the vicin- 
ity of bifurcating orbits. The scaling and uniformization predictions based on the 
semiclassical theory describe the variations of global quantities like the average in- 
verse participation number of the states. We would like to stress that the properties 
of the fixed points (e.g. through local Lyapunov exponent as in Hellers scar theory) 
do not explain the localization since the effect of the bifurcation is seen over a larger 
scale in the parameter K than expected from the appearance of the stable island and 
the local escape time calculated from the Lyapunov exponent. Such a theory could 
also not explain the interference pattern seen in the inverse participation ratio. 
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